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I. OVERVIEW 

The problem of explaining quark confinement has been 
of foremost importance since the formulation of quantum 
chromodynamics (QCD). The principal manifestation of 
confinement is the linear growth of QCD potential be- 
tween color charges. This is known to be the property of 
the Wilson loop l]. However, it has been impossible so 
far to use this in an analytic ab initio calculation in 3 + 1 
dimensional QCD. We want to deal with this challenging 
problem by combining: 

a) the Erickson-Semenoff-Szabo-Zarembo (ESSZ) [H, Q 
formulation for Bethe-Salpeter type equation for Wilson 
loops, with 

b) Dyson-Schwinger equations (DSE) for the gluon and 
ghost propagator in Landau gauge [j, [![ . 

We solve DSE for gluons and ghosts in the pure glue 
two-point sector. Then we insert the resulting QCD cou- 
pling a = g'^/iw and the gluon propagator into ESSZ 
equation for a rectangular (non-supersymmetric) Wilson 
loop, and solve this integral equation, which yields the 
Wilson potential. 

The ESSZ type ladder (or rainbow) diagram sum- 
mation has long been a major tool for extracting non- 
perturbative information about dynamics of a gauge the- 
ory. However the strength of this method is more evi- 
dent in A/" = 4 supersymmetric Yang-Mils due to higher 
order vertex correction cancelation. In principle the 
use of ESSZ ladder summation in our context of non- 
supersymmetric QCD is doubtful, and we will make sev- 
eral efforts to establish this approach: we will study the 
vertex correction terms by comparing the leading order 
(LO) contribution to the next to LO (NLO) contribution 
of the three-gluon vertex, and we will consider conver- 
gence of the entire procedure by evaluating the string 
tension at different DSE scale fixing points. 



Within QCD, the DSE for propagators and vertex 
functions have been studied in great depth, for review 
see [1, Q and references therein. Relation of DSE to 
lattice results is discussed in @. An alternative related 
method of functional renormalization group has been dis- 
cussed in 0. The relevant results on three-point func- 
tions are seen in Q , and on quark propagator in Q , the 
question of confinement inherent alone in DSE are dis- 
cussed in [13, [m , the uniqueness of the infrared (IR) 
scaling of Green functions established and gluon prop- 
agator IR non-singularity strictly supported in [12) > 
the IR universality established in [l4|. 

Below in section HIl we describe the ESSZ equations in 
a pure Yang-Mills theory with an arbitrary propagator 
(form- factor) . In section Hill we present DSE and our so- 
lution, our results are in agreement with the standard 
state-of-the art calculations of ghost and gluon propaga- 
tors in Landau gauge. In section lTVl we evaluate the ESSZ 
truncated Wilson loop employing the DSE propagators 
from section IIIII and check the significance of the NLO 
vertex correction. In section |V] we discuss the reasons 
why confining potential is not observed either in pure- 
glue two-point sector of DSE, or ESSZ solely, yet it is 
seen in the combination thereof. 



II. ESSZ EQUATION 

The Wilson loop 



W{C) = (tr Pexpj^ A^{x)dx^'j) (1) 

offers information about the behaviour of quarks in the 
theory, and the quark-antiquark potential is 



y(L) = - hm ilniy(CT,L), 

I — >oo 1 



(2) 
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where Ct,l is a rectangular Wilson loop in the {x^,x^) 
plane, with T being loop temporal length, and L loop 
spatial length, T ^ L. 
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The Wilson loop Eq. ([T]) can be represented in terms of 
a perturbative expansion, which can be found e.g. in the 
review [l^ . A set of Feynman rules for Wilson loops can 
be found in [l^, which will be of use to us below. Per- 
turbative treatment of Wilson loops is not useful in the 
non-Abelian case, and especially in the present context, 
as it yields obviously wrong results for Yang-Mills theory, 
for which it predicts a Coulomb-type potential [l^. A 
large- A'^c partial summation of ladder diagrams has been 
proposed in 0, 0] and performed for a circular and a rect- 
angular loop in J\f — A supersymmetric model (SUSY). 
This method is adapted here to the case of a non-SUSY 
theory, for the case that the partial summation of per- 
turbation theory (PT) series for propagators has already 
been performed in terms of solving DSE. 

Consider a trapezoidal loop W{C) — T{Ti, T2; L) with 
long parallel temporal sides of lengths Ti,T2, separated 
by a spatial distance L. Then the requirement that 
adding a propagator to the summed expression does not 
change it leads to the following integral equation for the 
sum of all ladder diagrams: 



r(Ti,T2,L) = 1 



47r2 



Ti 



dti 



T2 



dt2 X 



(3) 



xT{ti,t2,L)D^,{{xi - X2f)x';x^2, 

dots denote derivatives in ti^2 respectively, where x^ = 
Xi{ti), X2 — 2^2(^1) paths running over the Wilson 
loop as functions of ti,t2. For a rectangular loop xi — 
(— L/2, ii, 0, 0),a;2 = (L/2, t2, 0, 0). Configuration space 
propagator is related to the momentum-space form-factor 
F{p'^), introduced in the next Section Hill bv: 



D^,{x^ 



(27r) 



■tpx 



p2 ; 



(4) 

For simplicity we write Df^^{x'^)x^X2 = D{x^). Bound- 
ary conditions imposed upon F are 



F(r,0;L) = r(0, T;L) = 1. 



(5) 



The potential is related to r{Ti,T2; L) in the following 
way 



V(L) = - lim -lnF(r,r:i) 



(6) 



Equation ^ is depicted symbolically in Fig. Obvi- 
ously, if we write down the first term for T{Ti,T2] L) in 
the expansion of the solution, we shall reproduce the 
perturbative result for the Wilson loop. 

The central filled square in Fig. ^ symbolizes an ir- 
reducible kernel, containing (potentially) all the possible 
loop corrections. A convenient way of solving Eq. ([3]) is 
to consider the equivalent differential equation: 



dtidt2 



47r2 



(7) 



Wilson loop 




1 + g N^. idt,dt, 



FIG. 1: Summation of ladder/rainbows for a Wilson loop. 



We now introduce the variables x — {ti — 12)/ L, y ~ 
{ti +t2)/L. With this Ansatz the separation of variables 
becomes possible, and using the form: 



F = y^-0n(a;)e 



2L 



we will be solving the Id-equation 



^^n{x) +U{x]L)^Jn{x) 



dx 

with the effective potential 



(8) 



i^n{x), (9) 



U{x;L) = 



5^c,2 

47r2 



L^D{L\\+x^)) 



(10) 



We are solely interested in the unique ground state so- 
lution of Eq. ([5]), since the Wilson quark-quark potential 
is 



V{L) 



(11) 



A degeneracy in solutions of Eq. © may arise and thus 
complicate the situation, however, we have never ob- 
served it in our numeric calculations shown below in sec- 
tion |TVl It is now evident, that in order to complete the 
Wilson potential evaluation we need the propagator D 
and the coupling a = g'^/in derived from DSE in order 
to be able to evaluate V — Qo{L)/L. 



III. DYSON-SCHWINGER EQUATIONS 

We now obtain the nonperturbative input to ESSZ 
equations, i.e. the Dyson-Schwinger improved gluon 
propagator and coupling a. The difference between DSE 
and the simple renormalization group (RG) improved 
quantity is in the IR and medium momentum ranges, 
their ultra violet (UV) behaviour being identical (up to 
1 loop at least). Our DSE procedure uses the technique 
described in [3, US] > the reader familiar with this may 
skip the current section where we demonstrate that the 
results of [l^, [lO] are independently reproduced by us. 

We employ in our work the Newton-method based nu- 
merical technique described in [2l|. We solve a system 
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for ghost and gluon propagators, corresponding to the 
representation seen in figure Fig. Here bulbs de- 

note dressing of the propagators, and transparent bulbs 
- dressing of vertices. 



-1 



P -1 



1 q 



Ghost 

P-1 P -1 

V\/^j\f\ = \AAAAA/\ 

Gluon ^0' 
FIG. 2: Diagrammatic representation of DSE. 

These equations can be written in the form: 



G(p2) G(/I2) 

1 1 



(12) 



(n(p2)-n(/22)) 



where vacuum polarization is 

U{p') = U'^{p')+U'\p'), (13) 

tf^(/) = N^g^J ^^M,{j>\q\r^)G{q^)G{r% 



and self-energy is 



np')^N^9' / K,{p\q\r^)G{q^)F{r^)0^ 



(14) 
(15) 



Here Jlg^c are subtraction points, p,c = 0, = M, M is the 
limit of the interval p^ € (0, p?) in the momentum space 
where we solve DSE, coupling is meant to be g'^{p,^). F 
is gluon propagator form-factor in Landau gauge, defined 
via relation 



p J P + *e 



(16) 



and the ghost propagator non-trivial behaviour is de- 
scribed by the form-factor G 



D^^\p)^^Gip'). 
-I- le 

Variable z is the logarithmic variable 



2; = In 



El 



(17) 



(18) 



and scale /z is yet to be defined upon solving Dyson- 
Schwinger equations from comparing the obtained cou- 
pling aDSE{z) to the known values of apoaip^) at point 
M: 

aDSEiHMyfi^)) = apDciM^)- (19) 

The coupling constant 5^/47r = a is expressed in terms 
of G, F solely as vertex is finite in Landau gauge 

(at one-loop level) 

aDSE(ln(p')) = aBSEia)F{p^)G\p''). (20) 

In our case, we shall use varying scale fixing point M 
so that we can prove that our results are independent of 
scale fixing point choice within the error margin of our 
procedure. 

The kernels Mq, Kq,Qq arc known in literature, but 
for self-containedness of the paper we show them here: 



Ka{x,y,e) = 
Mo{x,y,6) = 



sin-^(9) 



(21) 



12x {-2coa{d)^/xy + X + y)' 



{ysm^{9) [2cos(26') {'ox'^ + ilxy + &y'^) - 
— 12xcos{?>9),/xy + xycoa{46) — 48 cos(6')y/x?7(x -I- y) — 
-12y cos{39)^/mj + 3x'^ + 27xy + Sy^] } . 

(22) 

For convenience, variables x = p^, y = q^ are introduced; 
variable 9 is defined via {p — q)^ — x + y — 2Y^a?y cos 9. 

To solve Dyson-Schwinger equations we use the 
Ansatz [11, [H: 



exp j aiTi{z) ] , z £ (lne,ln/i^ 



F{z) = I 



J2\ 1 



G{z) = 



F{fl) {l+uj\og^] ,z>ln^^ 



Az^", z < Ine, 



exp I ^5iTi(z) j , z e (lne,ln/i^ 



(23) 



G((t) 1 +tjln^ ,z > ln/2 



Bz z < e. 



Here are Tschebyschev polynomials, 0^,6^ are un- 
known coefficients yet to be determined from the nu- 
merical solution, n is the number of polynomials used 
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(mostly n = 30 has been used here, ahowing precision 
of 10-1'' for coefficients), 6 = -9/44, 7 = -1 - 25, 
ijj — llA^r.a((7)/(127r). The IR scahng k is chosen to 
be the standard 122, l23l 



K = 0.59 



(24) 



for the case of Brown-Pennington truncation with Q = 
1 [iO] (for discussion of meaning of C, see [l^l), which is 
our case (Calready set to its number value everywhere). 
Following [5| , we employ renormalization constant Zi re- 
definition, so that no momentum dependence could pos- 
sibly enter it, that is 



(^(y)(l-a/5-2a) Q(^y^^(l-b/ S-2b) 

^1 = F(y)(i+°) F(?/)(i+&) 
Again, following Q we choose 

a — b ~ ^iS, 



(25) 



(26) 



which minimizes its momentum dependence. Renormal- 
ization constant Zi refers to the piece with ghost loop in 
vacuum polarization. The equations are solved by using 
Newton's method, very clearly described for this particu- 
lar application by Bloch [2l| . The results of the solution 
are propagator formfactors F, G, shown in Fig. ([3]) , the 
IR behaviour of the propagators corresponds to the stan- 
dard ghost enhancement and gluon suppression. 

log,oG(p') 



In(pV) 




logioF(P ) 



FIG. 3: Ghost (dashed line) and gluon (solid line) propagator 
formfactors obtained in DSE in Landau gauge. 



The coupling a obtained from DSE Eq. (|20p is shown 
in Fig. (HI). We compare it to the standard coupling from 
Particle Data Group [2^ , and note that the both coincide 
very well in the UV. We also note here that the IR fixed 
point seen in the Figure is 

a(0) w 2.9 

for Nc = 3, which is consistent with the up-to-date 
Dyson-Schwinger results reported by other groups [1, Q . 




In(pV) 



FIG. 4: Running coupling from Dyson-Schwinger equations: 
our solution of DSE (solid line), and standard PDG coupling 
(dashed) . 



IV. SOLVING ESSZ EQUATION 

We have to find the lowest eigenvalue of a Schrodinger 
equation Eq. © 



"2dx2 



+ U{x;L) V(a;) = S^^ix) 



(27) 



where the auxiliary potential U{x) is related by an linear 
integral transform to gluon form-factor as 



U{x) = -- 



2'KaNr. 



1 



du 



(l + a;2)2 (27r)2 
(uJi(m) - (1 - 2,x^)J2{u)) F ( In 



(28) 

where /i is defined at point M as given in (fT9|) , M varying 
from 1 to 10 GeV, u is dummy scalar dimensionless in- 
tegration variable. The coupling a, in the sense of DSE 
approach, is taken here at the scale of rather than 
bare. We solve the Schrodinger equation with shooting 
method and find its ground state. Special care is taken 
to make sure this state is not degenerate. As a result we 



get the QCD potential V{L) 



The potential is 



defined up to additive constant, so we shift it to provide 
convenient comparison to existent results. It is shown 
in Fig. (O below, and is compared with lattice results by 
Gubarev et al. [2^ and Necco [23|. Linear IR behaviour 
of the potential can be clearly seen from Figure. We fit 
the potential by the standard expression 



V{L) 



4 oq 
3 L 



- Co 



(tL. 



(29) 



Dependence on string tension a on the scale fixing point 
choice is shown in Fig. ©. We see that the variance of a 
does not exceed that of different lattice results, shown in 
the table (jT|. The error we quote arises from an average 
of results obtained at different scale fixing points. This 
yields ao = 0.24 and a = 1.07 ± 0.1. 
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TABLE I: Comparison of string tension from different sources 



V[GeV] 



Author 


Year 


cr,GeV/fm 


Bali et al. [28] 


2000 


f.27 


Necco [27,1 


2003 


f.f9 


Gubarev et al. 


26] 


2007 


0.978 


Weise et al. \29^ 




2009 


f.07 


Present work 


2009 


f.07±0.i 





FIG. 7: Terms in ESSZ equations generating two-loop correc- 
tion to the auxiliary potential U'^^^x). 



as defined in llal we obtain: 



FIG. 5: Potential as function of distance, solid line - our 
result, dashed line - result by Necco, 2003 [131, dotted line - 
result by Gubarev et al. 2007 Hi]. 



The key result, the linear confining potential comes 
as a surprise. It invites the question, how great are the 
corrections coming from three-point vertex? One actu- 
ally shouldn't have thought that QCD can be described 
with ESSZ partial summation structure. Considering the 
vertex, the auxiliary potential is then modified: 



U{x) ^ U^^^x) +4TTaNcU^^\x), 



(30) 



where C/'^'(a;) comes, in the leading 1/Nc order, from the 
Wilson loop diagram shown in Fig. ([7]). 

Calculating the diagram in Landau gauge with rules 



n[GeV/fm] 
hl5 



1.05 



1.00 




M[GeV] 



FIG. 6: Dependence of string tension a on scale fixing point 
M. 



{uiU2){u^y) 
with 



1 



{y-xiY{y-x2)'^{y-x:iY 
1 1 



{y-xiY {y-X2)' 



cyclic permut. 

(31) 



f XI = (-L/2,ti,0,0) 
X2 = (L/2,t2,0,0) 

Xi = {-L/2 + Lts, h + h{t2 - h), 0, 0) 
ui = (0,1,0,0) 
U2 = (0,1,0,0) 
L M3 = (i,t2-ii,0,0) 



(32) 



Numerical evaluation of this integral shows that within 
the whole range of values of ii,t2 with which we work, 
U^'^\ti,t2) = C/(2)(ti - t2) = U^'^\x). This makes sep- 
aration of variables still possible and provides an extra 
test for reasonability of our model. Numerical values 

of U^^^{x) are such that U'^'^\x) /U'^^\x) w 10-^, which 
makes its contribution to auxiliary potential ground state 
negligible. This allows us to justify validity of ESSZ equa- 
tion application in the non-SUSY case: vertex correction 
is present but numerically suppressed. 



V. DISCUSSION AND CONCLUSION 

ESSZ approach to SUSY Wilson loops has worked very 
well in [3, 0]. The reason for that is absence of NLO 
corrections in the maximally supersymmetric theory. At 
small coupling their result has restored the perturbatively 
known IR singularity structure. Moreover, the calcu- 
lation originally performed in the small coupling limit, 
could be continued into large coupling limit. At large 
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coupling the solution to ESSZ equation reproduces al- 
most exactly the gravity dual result [H, [s^, HI] (up to 
an overall numerical factor very close to unity). Actu- 
ally this result, though in a different theory, has been a 
guide for our QCD treatment: as we are dealing with the 
IR strongly coupled theory, we are certainly out of order 
of applicability of any perturbative treatment, and even 
summation of diagrams would be suspicious. 

The reason why the ESSZ equation has never been 
applied to non-SUSY contents is obvious. It is clear 
from 2, 3] that when a perturbative propagator input 
is being used only a non-confining Wilson loop, with a 
Coulomb-type potential may be obtained. This follows 
from the fact that dependence on Wilson loop spatial size 
L may be scaled out of the ESSZ equation, so that any 
potentials one gets from it are Coulombic, varying from 
each other by coupling rather than distance dependence. 
Thus such a result would have been a priori useless in un- 
derstanding anything about strong coupling IR regime of 
gauge theory, where confinement governs the dynamics. 
This maybe the reason why summation a la ESSZ has 
not before been employed in pure Yang-Mills theory. 

A direct perturbative calculation of Wilson loop with 
a Dyson-Schwinger propagator yields no confinement 
whatsoever. Only quark-gluon vertex functions P.[Tol.[Tll| 
coming from DSE can render something looking like con- 
finement, which is then related to the singular behaviour 
of the quark-gluon vertex in the IR. 

In our opinion, this could not be one of the possi- 
ble ways to approach within the DSE the confinement 
problem, since it requires quark coupling to be singular, 
gluon one regular, which constitutes a severe violation of 
Ward identities. Rather than to involve three-particle 
functions, we apply ESSZ summation with DSE solu- 
tions which are possessing intrinsic scale, distance L is 
no more possible to scale out of ESSZ equations. Thus 
the resulting Wilson potential is no more necessarily be- 
ing Coulombic. 

A description of a single Wilson loop, from which one 
can obtain the QCD potential and provide a criterion of 
confinement, has not been done so far in terms of the 
two-point sector of DSE hierarchy. Thus our work closes 
an essential gap in the literature. The main reason for 
this gap was the theorem by West [l^l , stating that con- 
finement is provided by a very IR-singular propagator 
D{q^) ~ 1 /(?''. We know however that gluon propagator 
is regular in the IR in the DSE approach. 

Our work is based on combined analysis of Green func- 
tions and Wilson loops, allowing thus a study of the spa- 
tial QCD potential. This distinguishes our approach from 
several earlier papers where gluon non-propagation was 
considered instead of confinement and related to the an- 
alytic properties of Green functions, in particular, to the 
IR scaling k, Eq. ([M)) . These other works use the word 
"confinement" as in the original paper [ij when they 
mean to say of "non-propagation" . Known are the so- 
called Kugo-Ojima criterion for colour non-propagation 
K > [30| , Zwanziger criterion of ghost non-propagation 



K > and gluon non- propagation k > 1/2 [3l|. A claim 
has been made |32] for k > 1/4 to be quark confinement 
criterion by analysis of the Polyakov loop and effective 
QCD action in an external field. All these results are 
about gluon non-propagation rather than the properties 
of a colour charge-colour charge confining interaction. 

Returning to the discussion of our results we note that 
the reliability thereof may be questioned in what con- 
cerns the DSE input. The first issue is the truncation 
of the DSE system we solve to only two-point func- 
tions. The truncation is justified by ghost-gluon ver- 
tex not acquiring acquire one-loop corrections in Landau 
gauge. It has been proved that the three-point gluon 
and quark-gluon functions don't change ghost dominance 
property [8|, even though they are important for bound 
states |37| . In this sense, vertex functions are unimpor- 
tant for our particular context. 

Another question is whether Green functions obtained 
from DSE are physically relevant within the Wilson loop 
context we are discussing. We note that it is mostly 
medium-energy range that provides the important con- 
tribution into the auxiliary potential U{x; L), rather than 
the perhaps more model dependent IR piece. The Wilson 
loop thus depends on medium energy range values of the 
propagators where the DSE behaviour is the same as in 
lattice. There are unresolved questions regarding com- 
parison of IR scaling |38l] within lattice and DSE. These 
issues have yet to be understood and resolved, although 
they do not affect our results materially. 

The observables a, a we compute are in principle gauge 
invariant. Our results are obtained in Landau gauge, 
which, as noted, is a convenient choice. It should be 
possible to check gauge-invariance explicitly at one loop 
level, we however do not do that here, since this tran- 
scends the scope of the present paper. We think that the 
possibility for the observable we consider to be gauge in- 
variant at one-loop level comes from the fact that several 
gauge-dependent objects are combined. 

We speculate here d propos that a nonperturbative 
summation a la ESSZ could improve significantly the 
properties of a correlator of gluon strengths with Wil- 
son lines 

J'ix) = {tYF^,{x)U{C)F''''{0)U+{C)) (33) 
U{C) being a phase factor 

U{C) ^Pexpl^ig J A^dx^j, (34) 

which differs from Wilson loop since the path is connect- 
ing the arguments in Eq. ([55)) i.e. points x and 0. Eq. ([55]) 
hade recently been of great interest [3§| as it represents 
an important vacuum property. As far as we know, a 
Bethe-Salpeter equation for this kind of correlator has 
not been developed yet. We attempted to evaluate it 
perturbatively [4Q|. The present effort arose from this 
earlier one but should have actually anteceded it, for then 
a framework for ESSZ summation may have been closer 
or even in hand. 
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A hypothesis should be considered that using a rel- 
evant component of the non-perturbative input from 
Dyson-Schwinger equations one may be able to obtain 
a self-consistent picture of the QCD vacuum with all 
higher correlation functions, colour confinement and con- 
densates, which is supported at the simplest LO level by 
the presented calculation. 

To conclude, combining Dyson-Schwinger summation 
for gluon and ghost propagators with the Ericson- 
Semenoff-Szabo-Zarembo summation (truncation) for 
Wilson loop, we have obtained the string tension and 
have further demonstrated that its value is nearly not 
dependent on the selection of the DSE scale fixing point, 
thus establishing the internal consistency of this novel 
description of confinement. The string tension deter- 
mined by our method for the pure SU(3) gauge theory is 
a = l.OTibO.l. The UV Couloumb behaviour is governed 
by ao ~ 0.24. 

One can actually be quite amazed that our method has 
worked so well in QCD, without supersymmetry, thus 
with vertices non-compensated. One can speculate that 
the two truncated summations are complementary, ESSZ 



taking care of ladders and DSE taking care of rainbows 
in the vertices. Among interesting further steps in the 
development of this framework we recognize the formu- 
lation and evaluation of a similar ESSZ equation or a 
correlator of two giuons, having in mind its application 
to gluon non-local condensate Eg. ([55]) . Another, per- 
haps more challenging further development could be to 
solve ESSZ and DSE jointly, without the separation into 
partial systems. 
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